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Abstract: We present the simplest model for classical transitions in flux vacua. A 
complex fleld with a spontaneously broken U{1) symmetry is embedded in M2 x Si. 
We numerically construct different winding number vacua, the vortices interpolating 
between them, and simulate the collisions of these vortices. We show that classical 
transitions are generic at large boosts, independent of whether or not vortices miss 
each other in the compact Si. 
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1 Introduction 

Classical transitions [1-6] — first order phase transitions facilitated by domain wall collisions 
have recently been studied in various situations. One of the original motivations is the 
simple intuition in [7] about fiux vacua. As shown in Fig. 1, different fiux vacua have 
diflFerent numbers of field lines threading through the compact extra dimensions. The 
boundary between them must contain charges — objects at which field lines can origi- 
nate or terminate. If one integrates out the details in the compact dimensions, these 
boundaries are just domain walls. However, as localized objects in the extra dimen- 
sions, a collision between domain walls does not imply a collision of these charges. 
When the charges reside in different places in the extra dimensions, they may simply 
miss each other. The natural consequence of this "missing" is leaving behind a region 
with less fiux — a transition into a lower vacuum. 

There are of course complications to such a simple idea. For one thing, both charges 
and fiux lines may have strong back reactions on the geometry [8]. Even with a fixed 
background geometry, the charges will have long range interactions with each other. It 
is not certain that a clean picture as in Fig. 1 is realistic. 

On the other hand, it is known that solitonic objects tend to pass through each 
other if boosted to large velocities, even if they collide head-on. This means that 
classical transitions can take place when vortices collide. In this paper we present high- 
resolution simulations that further clarify the dynamics of classical transitions between 
fiux vacua in this model. 
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Figure 1. (left) A conformal diagram showing collision of two domain walls (red and green 
lines) where the compact extra dimensions is suppressed. Initially they separate regions con- 
taining 3 fluxes and 2 fluxes as labeled. After the two walls pass through each other, they end 
up separating regions containing 1 and 2 fluxes, (right) We chose two time slices, taken at 
the two times deflned by the dashed lines on the conformal diagram, in which we include the 
compact extra dimension. The blue flux lines can only originate or terminate at charges. The 
red and green domain walls are actually point-like charges, represented by circles and squares 
respectively. The charges start in the state deflned by the lower tube, representing an early 
time, then follow the arrow of the corresponding color and evolve to the state deflned in the 
top tube. It is simply Gauss' law that demands the middle region being left with one flux, 
resulting in a different compactifled vacuum. 



2 Simplest Model of Flux Vacua 



We adopt the simplest model in [7]. Consider a complex scalar field (j) with a sponta- 
neously broken U{1) symmetry, 



We embed this theory in a (1+1+1) dimensional spacetime, M2 x S^i, 



S = I dt dx dy 



^(0) , 



(2.1) 



(2.2) 



where < ^ < L is the coordinate for the compact extra dimension. This leads to a 
discrete family of vacua, 



exp 



27my 



(2.3) 



Here the integer n < is the winding number of the vacuum. This also defines the 



energy density spectrum, 



E„ 



2 ' 



L2 
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The phase < ^ < 27r is an exact symmetry and means nothing for a single vacuum 
state. However, when we put two vacua together, 



(j){x,y) = (j)n^eL^ for x^-oc, 

y) = (p(n+i),eR^ for x ^ oc, (2.5) 

the relative phase 

Ae = {eL-en), (2.6) 

plays an important role. The interpolation of minimum energy between these two 
asymptotic forms must contain a single vortex — a localized region of high energy 
density. Since their phases agree at 

this phase remains constant through out the interpolation. Therefore, at the opposite 
side, 

y = L-^. (2.8) 

the phase must go through the most dramatic change. That is the center of the vortex. 

This vortex will be attracted to the (n + l)th vacuum, and eventually relax the 
entire region into the lower vacuum. If we integrate out the compact dimension we 
describe this is a domain wall in the x — t spacetime that moves to engulf the region of 
higher vacuum energy. To produce collisions of domain walls, we can setup the initial 
condition with three regions, 

left, 

middle, (2.9) 
right. 



A6> 
2 



(n+l),0 



A6> 
' 2 



From the x — t point of view, it is a shrinking false vacuum region with two domain 
walls that slam into each other. However when 7^ 0, it seems like the two vortices 
will not collide. This is exactly what the pictorial understanding we gain from Fig. 1 
leads us to believe; hence, we might expect classical transitions to occur whenever the 
localized vortices do not collide. 

However, we should be careful before exercising such intuition. Saying that the 
vortices miss each other implies that they are much smaller than L. It is well known 
that in a non-compact 2D space, the vortex of this global U{1) does not have a finite 
size. That fact of course remains after we compact ify one of the spatial dimensions. 
Although the center of two vortices may have a large y separation, their profiles still 
overlap significantly during the collision. 
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On the other hand, overlapping profiles does not guarantee failure of classical tran- 
sitions. In[l, 2] there was no extra dimensions and the domain walls simply collide with 
each other. It was shown that when the collision is very energetic, the field profiles are 
so highly boosted that they cannot efficiently interact, therefore must remain intact. 

Given these two motivating reasons to believe classical transitions should occur in 
our model, one might expect that large separation, or highly boosted domain walls 
will lead to classical transitions. At the same time, there should be a low boost and/or 
small separation regieme where it fails. The result turns out to be more subtle, as we 
shall demonstrate later. 

3 Relaxing to the Wall 

The first challenge that we need to address is how to set up the initial conditions. 
The functional form for the field profile that interpolates between local minima of a 
potential [9-11] is well known for some potentials. In some cases these initial condi- 
tions are analytic; in others relaxation to an approximate solution is necessary before 
simulations can be run. In the model presented here, our vacuum states are defined 
by field configurations in the compact dimension and the interpolation between those 
states must be found numerically [12, 13] by relaxing the field to the instanton solution. 

The vacuum states of our model are defined by Eq. (2.3). Our next goal is to 
numerically setup the situation in Eq. (2.9) that will lead to domain wall collisions. 
For numerical reasons, even the uncompactified dimension x can only have a finite size. 
We will take that to be —Lx/2 < x < Lx/2 and enforce fixed boundary conditions at 
X = ±Lx/2. One should not confuse the finite box size Lx with the size of compact 
extra dimensions L in Eq. (2.3). We will maintain that Lx ^ L through out the paper. 

The initial field profile is a smooth interpolation version of Eq. (2.9) with n = 1. 



Fig. 2 shows the real and imaginary parts of Eq. (3.1) along with a plot of the energy 
density. 

Relaxation is achieved by imposing a period of high friction at the start of the 
simulation which damps out spurious radiation produced by the field as it relaxes to 
the instanton solution. The field values for the majority of the volume of the simulation 
should be very close to the vacuum states, Eq. (2.3), with localized, although not 
infinitely so, vortices along the domain walls that separate the different states. The y 
separation of the two vortices should follows from our setup of A9. This is indeed what 



4>it = 0) 




cos' [t^j 4>2fi + sin^ i^f-x 
cos^ ( -^x ] 02,0 + sin^ ( -^x 




(3.1) 
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y y y 

Figure 2. (left) The real, (f)R, and (middle) imaginary, 0^^, parts of the initial field profile, 
(right) The energy density of the initial field profile. 



we see in Fig. 3. These define the initial conditions for the next step: acceleration and 
collision. 




y y y 



Figure 3. (left) The real, and (middle) imaginary, 0^^, parts of the field profile after 
relaxation, and also the energy density. We can see from the energy density plot (right) that 
it is highly concentrated at two points, which are basically the center of vortices. We can also 
see the three vacuum regions in which the energy density is almost constant, and the small 
difference in energy between the middle region and the ends. 

Note that, in principle, the motion we hope to see is already happening during 
this period of high friction. We are just exploiting a separation of scales. The even- 
tual accelerations of vortices happens in a longer time scale, comparing to the energy 
damping from Fig. 2 to Fig. 3. In practice, we execute both relaxation to the initial 
condition and the evolution afterward in one numerical process. The simulation begins 
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at t = 0, and reaches Fig. 3 at t = t^, at which point the friction is turned off (or 
greatly reduced). 

We shaU also introduce the "Einstein Frame" point of view, in which the details in 
the compact dimension y are integrated out. Fig. 4 shows the energy density integrated 
over the compact dimension. At time tg we see that the high peaks in Fig. 3 are actually 
pretty low energy barriers between the two vacua. That is because the n = 2 state in 
the middle is the highest vacuum in this parameter choice. We also present the curves 
at later times when the walls pick up some velocity. We can see the integrated energy 
density of the wall acquires a Lorentz boost and grows in magnitude. This confirms 
that there is indeed some "object" being accelerated during the simulation. 



Figure 4. The integrated energy density as a function of x. The thick-short dashed curve is 
at t = tg. The solid and dashed curves are at subsequent times. 

4 Collisions 

After turning off friction at t^, the field continues to evolve according to the equation 
of motion one can derive from the action in Eq. (2.2). In practice we still keep a small 
non-zero friction to control numerical noises. As long as the resulting energy loss is 
tiny compared to the energy scale involve in the simulation, we claim that our result 
is a good approximation to the exact equation of motion. 

Before the collision, as shown in Fig. 4, the vacuum regions and domain walls are 
quite sharply recognizable just from the energy plot. However during the collision the 
non-linear dynamics of the field make energy plots less informative. We need another 
quantity to keep track of the dynamics. Recall that different vacuum states, Eq. (2.3) 
has different winding numbers. So we can use the general definition of winding number. 




c 



X 




(4.1) 
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to recognize different states. Take the closed contour C in Eq. (4.1) be a path along 
the compact dimension every x position has a well-defined value of w unless a vortex 
is exactly there. Keeping track of the value of it; at x = 0, and evaluating it at the end 
of the simluation, allows us to tell whether a classical transition has happened. 

We are interested in whether the existence of a classical transition depends on the 
collision energy, and the distances between vortices in the y dimension. The later can 
be controlled by A6^ in the initial conditions, Eq. (3.1), as explained in Sec. 3. The 
collision energy is controlled by the wall separation at t^, since they basically come 
from the potential energy in the middle region. We adjust this by changing without 
changing the definition of the initial profile in Eq. (3.1). The resulting wall separation 
can be evaluated numerically, and it is accurately proportional to Lx. Scanning through 

and wall separation, we present the results in Fig. 5. 

There are several features to point out in Fig. 5. First, as the collision energy 
increases, eventually classical transition happens for all Actually for any fixed 
increasing collision energy always helps^. Most interestingly, for a fixed wall separation, 
increasing A0 does not really help. Actually it is the opposite. All the head-on collisions 
lead to classical transitions. It starts to fail at some large A0, and eventually recovers 
at the maximum A9 = n. 

LO 
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Figure 5. Contour plot showing the final state of the simulation at x = 0. White dots (on 
grey) represent simulations in which the stable transition to n = occurred at x = 0. Black 
dots (on white) represent regions that did not transition to the n = state, and the final state 
at X = was n = 1. 

We select a few examples to show the winding number plots during the simulations. 

^With one exception when the wall separation is close to L. It is unclear what happens there. 
When the wall separation is comparable to the compact dimension size, our usual intuition might be 
misleading. For example it might not be true that larger separation means more collision energy there. 
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To understand the effect of changing the coUision energy, we keep = 97r/16 and 
show the case of waU separation x ^ 1.5L in Fig. 6, and ^ 1.9L in Fig. 7. That is the 
difference between making the transition or not. In Fig. 6, w{0) does dip down to for 
a moment, but we did not get two stable vortices flying apart to maintain and expand 
that region. At a certain moment in Fig. 7 we actually have 4 vortices separating the 
entire simulation range into 5 regions. After that the middle n = 1 region collapses 
and gives us a successful transition. 

r = rv + 1 1 .7 T = rv + 13.2 T = + 273 



E E 
Z Z 
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Figure 6. w{x) just before the collision (left), during the collision (middle) and the final state 
after the collision (right). This is for wall separation ^ 1.5L and = 97r/16. 
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Figure 7. The similar set of three snapshots of w{x)^ for wall separation ^ 1.9L and AO = 
97r/16. 

Next we keep the wall separation at ^ 1.9L and present = 0, ll7r/16 and tt. As 
shown in Fig. 8, the transition process at and tt are similar to each other, and much 
cleaner than Fig. 7 where AO = 97r/16. The non-transition result at ll7r/16 looks quite 
similar to Fig. 6. 



5 Discussion 



We have demonstrated numerical simulations of vortex collisions in the background with 
compact dimensions. Since vortices are the natural boundaries for vacua with different 
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Figure 8. The similar set of three snapshots of w{x)^ for wah separation ^ 1.9L and = 
(top), ll7r/16 (middle) and tt (bottom). 

winding number, we see the expected phenomena of classical transitions. Our result 
agrees with the understanding from simpler models [1, 2] that more energetic collisions 
help to facilitate the transitions. However, as explained in [2] we should again advocate 
caution against a naive energetic argument that larger collision energy helps to make 
two non-interacting vortices. Since before the collision we have two fast approaching 
vortices, there is always enough energy to make two out-going ones. Also as shown in 
Fig. 7, providing more initial energy to the collision may just create complicated final 
states, like more vortices or sourcing radiative modes, instead of directly helping the 
transition. So one should think beyond the simple energetic argument. 

The more intriguing result is that against the intuition of [7] , larger vortex separa- 
tion in the compact dimension actually reduces the chance of transit. The near head-on 
collisions {A9 close to 0) almost always lead to classical transitions, and failures start 
to appear as A9 increases. So vortices "missing" each other actually does not always 
help. More precisely, we cannot define "missing", since the vortices do not have finite 



sizes smaller than L, and there are always "side-ways" interactions. The interaction 
basically transfers x momentum to y momentum. With less x momentum, the vortices 
fail to fly apart. 

This also explains why at = tt we always get transitions, and they look very 
similar to the cases where = 0. In both cases, the Z2 symmetry forbids the vortices 
from acquiring any y momentum. So they tend to exit the collision without having 
actually collided. Although tt happens to be the maximum separation of vortices in 
this model, we do not think the success of transitions should be attributed to this 
separation. 

Finally, in the intermediate values of we observed interesting behavior. For 
example the production of two extra vortices, as shown in Fig. 7, is the most unexpected. 
The short range interaction between vortices is an unexplored subject and we look 
forward to further developments: including more complicated compactifled geometries. 
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